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A PERTURBATION  SOLUTION  OF  ROTOR  FLAPPING  STABILITY 


Wayne  Johnson* 

U.S  'rmy  Air  Mobility  RID  Laboratory 
Moffett  Field,  California  94035 


Abe  tract 

The  stability  of  the  flapp'.ig  motion  of  a single  blade 
of  a helicopter  rotor  is  examined  using  the  techniques  of 
perturbation  theory.  The  equation  of  motion  studied  is 
linear,  with  periodic  aerodynamic  coefficients  due  to  the 
forward  speed  of  the  rotor.  Blade  pitch  feedback  pro- 
portional to  both  flapping  displacement  and  flapping 
rate  is  included.  Four  cases  are  considered:  small  and 
large  advance  ratio,  and  small  and  large  Lock  number. 
The  solution  of  this  problem  demonstrates  the  informa- 
tion which  may  be  obtained  using  perturbation  techniques. 
This  paper  discusses  mainly  the  results  of  the  analysis 
rather  than  the  mathematics  of  their  derivation,  with 
primary  emphasis  on  the  eigenvalues  (l.e. , root  loci) 
as  Indicators  of  the  system  stability  and  response.  The 
small  advance  ratio  results  in  particular  are  very  useful, 
being  valid  out  to  an  advance  ratio  of  about  0.  5.  it  is 
concluded  that  perturbation  theory  is  a powerful  mail  3- 
matical  technique  which  should  prove  very  useful  in 
analyzing  seme  of  the  problems  of  helicopter  dynamics. 

Nomenclature 
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MP 

Mp 

Me 

LHP 

RHP 

LHS 

RHS 

Re 

lm 

P 
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x 


flap  proportional  feedback  gain 

flap  rate  feedback  gain 

aerodynamic  moment  due  to  flapping 
displacement 

aerodynamic  moment  due  to  flapolng  rate 

aerodynamic  moment  due  to  blade  pitch 

left-hand  plane 

right-hand  plane 

left-hand  side 

right-hand  side 

r'ml  part  of  a complex  number 

Imaginary  part  of  a complex  number 

flap  motion  degree  of  freedom 

blade  Lock  number 

eigenvalue  or  root  of  the  system 

rotor  advance  ratio  (forward  speed  divided 
by  rotor  tip  speed) 
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v rotating  natural  frequency  of  flap  motion 

(centrifugal  and  structural  stiffening),  non- 
dimenslonallzed  with  rotor  rotational  speed 

$ rotor  azimuth  angle,  measured  from 

downstream 

Of  ) "the  order  of' 

( ) conjugate  of  a complex  number 

Introduction 


This  paper  considers  the  application  of  perturbation 
techniques  to  helicopter  rotor  dynamics.  Perturbation 
theory  has  been  well  developed  In  recent  years,  but  has 
not  found  much  application  to  rotary  wing  problems. 
Classically  helicopter  engineering  has  made  us.  of  the 
same  perturbation  theories  that  fixed  wing  engiroerirg 
has,  for  example  lifting  line  theory  and  engineering  beam 
theory  (both  require  a large  blade  aspect  ratio).  Another 
classical  example  is  actuator  disk  theory  (a  large  number 
of  blades  is  required).  These  theories  were  developed  on 
an  intuitive  basis  however,  and  the  more  rigorous  mathe- 
matical techniques  of  perturbation  theory  have  not  yet 
found  widespread  use  for  rotary  wings.  The  classical 
applications  are  largely  for  aerodynamic  problems;  the 
mathematics  of  those  problems  can  be  very  complicated 
however  because  the  equations  involved  are  highly  non- 
linear partial  differential  equations.  The  treatment  of 
dynamic  problems  can  be  more  tractable  since  only 
ordinary  differential  equations  are  involved  Problems 
with  constant  coefficient  linear  differential  equations  can 
be  solved  exactly  with  well  established  methods,  so  for 
these  problems  the  extra  effort  of  perturbation  theory 
may  not  be  justified.  On  the  other  hand  for  problems 
with  time  varying  or  nonlinear  differential  equations  the 
only  solution  procedure  generally  applicable  is  the 
numerical  integration  of  the  equations  of  motion.  How- 
ever, purely  numerical  solutions  are  not  entirely  satis- 
factory for  obtaining  an  understanding  of  the  physical 
character  of  the  system,  or  for  formulating  general 
design  rules.  Furthermore,  an  analytic  solution  for  the 
general  case  would  be  difficult  to  obtain  (if  possible  at 
all)  and  would  be  so  complex  as  to  be  hardly  better  than 
the  numerical  solution.  The  only  systems  that  can  be 
practicably  handled  analytically  are  those  involving 
linear  constant  coefficient  differential  equations.  Per- 
turbation techniques  are  available  which  are  methods  to 
study  time  varying  or  nonlinear  systems  such  that  at 
each  step  in  the  analysis  only  linear  constant  coefficient 
equations  must  be  handled.  Time  varying  or  nonlinear 
differential  equations  are  characteristic  features  of 
helicopter  dynamics  and  aerodynamics,  primarily  due 
to  the  rotation  of  the  wing.  Thus  the  possibilities  for 
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the  use  of  perturbation  theory  In  rotary  wing  problems 
are  very  extensive. 

This  paper  considers  the  stability  of  the  flapping 
motion  of  a single  blade  of  a helicopter  rotor.  This  is 
a single  degree  of  freedom,  second-order  system,  with 
analytic  aerodynamic  coefficients.  The  governing  equa- 
tion is  Linear  with  time  varying  coefficients;  it  is  given 
below. 


p - v\  = y j (Mp  - Mg)  P + (Mp  - KpM^)  p 


(1) 


where  the  aerodynamic  coefficients  are 
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aiv!  the  three  regions  are  defined  by 


region  (i);  0 < p sin  t < P 

region  (ii):  - 1 < n sin  t < 0 

region  (iii):  - p < p sin  <p  < - 1 

This  is  the  homogeneous  equation  for  small  perturbations 
of  the  flapping  motion  of  the  blade  about  an  equilibrium 
state;  the  derivation  of  this  equation  may  be  found  in  the 
literature. W P is  the  degree  of  freedom  representing 
the  blade  flapping  motion  perturbation.  The  equation  is 
nondimensionalized  with  the  rotor  speed,  so  the  time 


variable  is  the  azimuth  angle  # . v is  the  rotating 
natural  frequency  (nondimensionalized  with  the  rotor 
speed)  of  the  flapping  motion,  which  may  be  greater 
than  1.0  for  flapping  hinge  offset  or  cantilever  root 
restraint  of  the  blade,  y is  the  Lock  number,  defined 
by  y = pacR^/Ijj  (p  is  the  air  density,  a the  two- 
dimensional  lift  curve  slope,  c the  blade  chord,  and  R 
the  blade  radius);  ijj  is  the  equivalent  mass  of  the 
flapping  motion,  given  by  the  integral  over  the  span  of 
the  square  of  the  mode  shape  of  the  flapping  motion 
weighted  by  the  mass  per  unit  length  of  the  blade;  for 
the  rigid  flapping  motion  of  an  articulated  blade,  the 
mode  shape  is  proportional  to  the  radial  distance  from 
the  hinge,  and  so  ljj  is  Just  the  moment  of  inertia  of 
the  blade  about  the  flapping  hinge.  Kp  is  the  flap 
proportional  feedback  gain,  better  known  as  tan  ; Kr 
is  the  flap  proportional  feedback  gain.  A feedback  law 
A 9 = -Kp/3  - Kp/!  has  been  used  (A9  is  the  blade 
pitch  change  due  to  flapping  feedback  control),  p is  .he 
rotor  advance  ratio  (forward  velocity  divided  by  rotor 
tip  speed).  The  coefficients  M^,  , and  are  the 

aerodynamic  forces  on  the  blade,  hence  their  multiplica- 
tion by  y and  their  dependence  on  p . The  three 
regions  for  the  coefficients  reflect  the  influence  of  the 
reverse  flow  region  of  the  rotor  disk.  In  region  (i) 
there  is  normal  flow  over  the  entire  blade  span;  in 
region  (iii)  reverse  flow  over  the  entire  span;  and  in 
region  (ii)  normal  fiow  outboard  of  r = - p sin  ip  and 
reverse  flow  inboard.  Region  (iii)  is  encountered  only 
if  p > 1 . The  aerodynamic  coefficients  were  ob- 
tained using  a rigid  blade  motion,  and  should  properly 
be  changed  some  to  handle  a blade  with  cantilever  root 
restraint.  However,  the  major  effects  of  a cantilever 
root  on  the  dynamics  of  the  system  are  due  to  changes  in 
v and  y (both  are  increased,  v to  1.15  say  and  V 
to  about  5/3  the  Lock  number  based  on  the  rigid  mode 
inertia).  Since  these  are  free  parameters  in  the 
analysis  this  formulation  of  the  problem  should  give 
reasonable  results  for  all  rotors. 


This  equation  has  been  studied  numerically  in 
recent  literature,  primarily  in  the  context  of  Floquet 
theory,  (2)  which  must  be  used  because  the  aerodynamic 
coefficients  are  periodic  in  <p  if  p / 0.  However, 
purely  numerical  solutions  are  not  very  satisfactory  for 
obtaining  an  understanding  of  the  physical  characteristics 
of  the  system,  or  for  formulating  general  design  rules. 
An  analytic  solution  for  this  equation  in  its  general  form 
would  be  difficult  to  obtain  (if  possible  at  all)  and  would 
be  so  complex  as  to  be  hardly  better  than  the  numerical 
solution.  Time  varying  or  nonlinear  differential  equa- 
tions are  characteristic  features  of  helicopter  dynamics 
and  aerodynamics,  primarily  due  to  the  rotation  of  the 
wing.  The  oniy  solution  procedure  generally  applicable 
to  such  problems  is  the  numerical  integration  of  the 
equations  of  motion,  in  particular  cases,  when  there  is 
some  parameter  which  is  either  smail  or  large,  the 
techniques  of  perturbation  theory  may  be  used  to  obtain 
analytic  solutions.  Thus  it  is  appropriate  to  appiy  the 
techniques  to  the  study  of  rotor  flapping  stability. 
Perturbation  theory  has  been  well  developed  in  recent 
years  (e.g. , Refs.  3 and  4),  but  has  not  yet  found  wide- 
spread application  to  rotary  wing  problems.  This  paper 
considers  an  example  of  such  an  application. 
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All  tbe  detail*  of  the  mathematical  analyst*  are  pot 
appropriate  for  this  paper.  The  mathematical  techniques 
are  not  very  sophisticated  actually;  there  are  some 
tricks  to  be  learned,  but  the  standard  ones  work  for 
most  systems,  Including  this  one.  (*•*)  However  the 
mathematics  can  be  rather  long  and  cumbersome, 
especially  when  the  higher  order  solutions  are  sought. 
Furthermore  it  is  not  maintained  that  the  differential 
equation  studied  is  a true  model  of  rotor  dynamics; 
nonlinear  aerodynamics  and  coupling  with  pitch  and  lag 
motions  are  certainly  very  important.  The  purpose  of 
this  papei  j not  to  describe  the  mathematical  details 
of  perturbation  techniques  or  to  present  a study  of  true 
flapping  dynamics.  Rather  It  is  intended  to  demon- 
strate what  information  can  be  o’jtained  by  the  perturba- 
tion techniques,  and  to  explore  the  methods  which  are 
most  useful  for  rotor  dynamics,  so  helicopter  engineers 
will  be  able  to  decide  whether  to  use  these  techniques 
with  more  complicated  or  more  realistic  systems. 
Therefore  only  an  outline  of  the  method  is  presented, 
with  primary  attention  to  the  results  obtained.  For 
more  details  of  the  analysis  the  reader  is  directed  to 
Ref.  5. 

For  the  flapping  equation  there  are  two  parameters 
which  may  be  used  for  perturbation  quantities;  the 
advance  ratio  p and  the  Lock  number  Y . Then  there 
are  four  cases  to  be  considered:  small  and  large  p , 
and  small  and  large  Y . The  dynamic  problem  con- 
sidered here  is  the  question  of  rotor  flapping  stability. 
The  stability  of  the  motion  is  determined  by  the  roots 
or  eigenvalues  of  the  system  (there  are  two  for  this 
second  order  equation)  and  so  most  of  the  results  wi’.' 
be  concerned  with  the  roots.  The  equation  considered  is 
linear;  perturbation  theory  is  used  because  the  aero- 
dynamic coefficients  are  time  varying  (specifically, 
periodic)  for  forward  flight,  i.e.  , when  p Is  greater 
than  zero.  Each  of  the  four  cases  for  this  problem  will 
be  considered  in  turn  below.  First  however,  brief  dis- 
cussions will  be  given  of  the  perturbation  techniques 
used  and  the  behavior  expected  of  the  eigenvalues  of  a 
periodic  system. 

Perturbation  Techniques 

Fundamental  to  perturbation  theory  is  the  existence 
of  some  parameter  which  is  either  very  small  or  very 
large;  for  the  moment  represent  the  small  parameter 
(or  the  Inverse  of  the  large  parameter)  by  c . In  the 
present  problem,  it  is  desired  to  find  the  roots  of  the 
system,  which  means  Investigating  a solution  which  is 
uniformly  valid  over  long  time  periods.  The  appro- 
priate perturbation  technique  is  the  method  of  multiple 
time  scales.  This  method  assumes  that  the  behavior 
of  the  system  may  be  investigated  over  several  time 
scales,  # = en#.  The  time  scales  #n  are  all 
assumed  to  be  of  the  same  order;  then  for  ^ = c # the 
actual  time  # must  of  order  c , I.e.  very  large 
comparea  to  the  basic  scale  #q  = # . The  derivative 
with  respect  to  # becomes  then 

d a a 2 _g_ 
d#  = ot0  + ea<pi  + € a*2  + • • • 


so  the  use  of  these  time  scales  amounts  to  an  expansion 
of  the  time  derivative  as  a series  in  c Next  the 
dependent  variable  is  expanded  as  a series  in  c 

(3  = p (a  . a , a ,..)+£  p (*  , a . . .)  «... 

o 1 V V *2  i *i  ’ 

where  the  terms  Pq,  Pj,  etc.  are  all  assumed  to  be  the 
same  order,  and  depend  on  all  the  time  scales.  In  addi- 
tion, all  free  parameters  (e.g. , v and  Y for  the  small 
It  case)  are  also  expanded  as  series  in  c . Those  series 
for  P , d/d#  and  the  free  parameters  are  substituted 
into  the  equation  of  motion,  and  all  terms  of  order  £ 
are  collected  and  separately  equated  to  zero,  giving  the 
o(en)  equation.  With  the  above  expansion  of  d/d  4> . 
these  equations  are  partial  differential  equations  now. 

The  requirement  that  all  PR  be  of  the  same  order  for 
the  long  time  scale  behavior  of  the  motion  is  crucial  to 
obtaining  the  solution  of  the  partial  differential  equations. 
It  leads,  for  certain  values  of  the  free  parameters,  to 
critical  regions  characterized  typically  by  a reduction 
of  the  stability  of  the  system.  The  method  of  multiple 
time  scales  is  used  for  the  small  p and  small  Y cases. 
Details  of  the  method,  and  its  use  for  these  two  cases, 
may  be  found  in  Refs.  3 and  5. 

Often  an  equation  of  motion  is  such  that  in  the  limit 
e = 0 the  order  of  the  differential  equation  is  reduced. 
The  large  Y and  large  p cases  arc  of  this  type;  for 
example  letting  Y — » results  in  a first  order  equation 
(just  the  RHS  of  Eq.  (1)),  which  is  independent  of  Y. 
Such  problems  are  called  boundary  layer  problems,  since 
they  are  characterized  by  narrow  regions  in  which  the 
solution  changes  greatly.  The  solution  outside  such 
regions  may  be  found  by  use  of  a substitution  of  the  form 
P = exp(/^pdi)  followed  by  an  expansion  of  p as  a 
series  in  e : 
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This  main  solution  is  not  valid  in  certain  narrow  transi- 
tion regions  or  boundary  layers.  A basic  part  of  this 
perturbation  technique  is  methods  to  obtain  solutions 
through  the  transition  regions,  so  that  it  is  possible  to 
match  one  main  solution  to  another  on  the  other  side  of 
the  transition  region,  or  to  boundary  conditions  at  the 
base  oi  the  boundary  layer.  The  matching  procedure  is 
so  fundamental  to  the  analysis  of  boundary- layer  types 
of  problems,  that  it  gives  the  technique  its  name:  the 
method  of  matched  asymptotic  expansions.  Details  of 
the  method  and  its  application  may  be  found  in  Refs.  3, 
4,  and  5.  Again,  for  stability  problems  the  behavior  of 
the  solution  over  long  time  periods  is  of  interest;  the 
method  of  matched  asymptotic  expansions  is  used  to 
construct  such  a solution,  from  which  the  eigenvalues  of 
the  system  may  be  found. 

The  methods  of  perturbation  theory  are  best  de- 
scribed by  example.  Since  the  mathematics  involved 
can  be  long,  espe  ally  for  higher  order  solution'',  a 
detailed  description  of  the  procedures  will  not  be  given 
here.  The  'jader  is  directed  to  Ref.  5 for  the  detailed 
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consideration  of  the  rotor  flapping  equation.  When  first 
learning  the  techniques,  the  reader  might  want  to  con- 
sider simpler  equations  as  examples,  and  these  may  be 
found  in  Refs . 3 and  4 . 

Periodic  Systems 

The  root  loci  of  a constant  coefficient  system  typi- 
cally exhibit  behavior  in  which  two  roots  start  as  com- 
plex conjugates,  met*  the  real  axis,  and  proceed  in 
opposite  directions  along  the  real  axis.  The  existence 
of  periodic  coefficients  in  the  differential  equations  of 
motion  generalizes  this  behavior  so  that  it  can  occur  at 
any  ImX  = n/rev  or  n-*l/2/rcv,  not  just  at  ImA  = 0. 
The  property  of  the  solution  that  allows  this  behavior  is 
the  fact  that  the  eigenvectors  are  themselves  periodic 
(instead  of  constant  as  for  a constant  coefficient  system). 
The  analysis  that  demonstrates  that  periodic  systems 
show  this  behavior  is  called  Floquet  theory.  Thus  the 
following  behavior  of  root  loci  is  characteristic  of 
periodic  systems.  If  the  parameter  being  varied,  for 
example  the  advance  ratio  p , is  such  that  at  p = 0 the 
system  is  not  periodic,  then  the  roots  start  out  as  com- 
plex conjugates . As  p increases,  the  periodicity  of  the 
system  increases,  and  the  roots  move  toward  n/rev 
(or  n * 1 /2/rev)  lines  (remaining  complex  conjugate 
pairs  though).  At  some  critical  p the  loci  reach 
ImA  = n/rev,  and  then  for  still  larger  p the  frequency 
remains  fixed  at  n/rev  while  the  real  part  of  one  root 
is  decreased  and  that  of  the  other  is  increased.  The 
root  beine  destabilized  usually  crosses  into  the  RHP 
for  some  p , indicating  the  system  has  become  unstable 
because  of  .he  Influence  of  the  periodic  coefficients. 

The  Small  p Case 


described  by  boundaries  in  p,  and  Y^  (or  p2  ■"d  *2  • 

etc.).  The  case  K = 0 will  be  considered  first. 

R 

Order  1 Results 

To  order  1 the  roots  are  the  hover  roots,  given  by 


and  its  conjugate.  The  variation  of  these  roots  with  V 
for  v = 1 and  several  Kp  is  shown  in  Fig.  1 (the 
p = 0 loci). 


Order  p Results 

To  order  p the  roots  for  most  pfl  and  Yfl  are 
stiil  the  hover  roots: 
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To  this  order  there  is  no  effect  of  the  periodic  coeffi- 
cients. The  exception  is  when  lmA0  1/2  (the 
frequency  of  the  hover  root  is  near  1/2  'rev)  i.e.  v0 
and  such  that 


The  small  parameter  in  this  case  is  the  advance 
ratio  p . Setting  p = 0 gives  the  hover  case;  the  time 
variation  of  the  aerodynamic  coefficients  is  removed  and 
the  roots  are  easily  obtained.  To  order  p the  reverse 
flow  region  may  be  ignored,  and  the  expressions  for 
Mp  , Mp  , and  Mg  in  region  (1)  may  be  used  for  all  41  • 
Tc  obtain  the  roots  it  is  necessary  to  find  a solution 
v-iiid  over  long  time  periods;  thus  the  appropriate  per- 
turbation technique  is  the  method  of  multiple  time 
scales.  This  method  Involves  considering  the  behavior 
of  P when  $ is  of  order  1,  of  order  P * , of  order 
p"2  , etc.  P is  expanded  as  a series  in  p : 


P 


nPj  ♦ p2P2 


and  the  method  yields  a solution  uniformly  valid  for  all 
ii  . p and  Y are  also  expanded  in  p : 


v = y 


m v. 


2 

M >o 


2 

v = •'o  + * "i  + M v2  + ’ ' ' 


There  is  then  a critical  region,  with  a boundary  described 
by  the  following  equation. 


This  equation  describes  a narrow  band,  of  width  O(p)  . 
about  IitiAq  = 1/2  (which  by  Eq.  (4)  may  be  considered 
a line  of  v as  a function  of  Y ) . CXitside  the  criMcal 
region  there  Is  an  0(p)  change  in  the  frequency  while 
the  real  part  of  the  root  is  unchanged  from  the  hover 
root.  Inside  the  critical  region  the  frequency  is  iixed 
at  l/2/rev  while  there  is  an  0 (p)  change  in  the  damping. 
This  behavior  is  illustrated  by  the  ;i  =0.1  loci  in  Fig.  1. 
The  maximum  stability  change  occurs  at  Pj  = Yj  0 
(i.e.  when  p and  Y arc  such  that  the  frequency  of  the 
hover  root  is  exactly  l/2/rev),  where  the  root  is 
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This  is  done  because  it  is  the  characteristic  of  systems 
with  periodic  coefficients  that  for  certain  values  of 
Pq  and  Yq  there  are  stability  degradation  regions 
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The  critical  region  is  a region  of  stability  degradation, 
not  of  instability;  the  maximum  reduction  (and  enhance- 
ment) in  the  damping  is 


(7) 

which  is  an  0(p)  small  reduction.  The  system  remains 
stable  because  of  the  large  hover  damping.  Fig.  1 shows 
typical  root  loci  for  varying  y , with  p = 0 and 
p = 0.1.  The  point  D on  the  hover  locus  (Kp  = 1)  is 
where  the  locus  crosses  ImX  = 1/2.  As  p increases, 
since  this  point  is  at  the  center  of  the  critical  region  it 
receives  the  maximum  stability  change,  and  so  is  pulled 
out  to  the  point  B.  In  terms  of  the  > locus,  as  V In- 
creases and  the  hover  locus  nears  ImA  = l/2  , the  root 
has  an  Oya>  change  in  the  frequency,  pulling  the  locus 
toward  ImA  = 1/2  . When  the  locus  crosses  into  the 
critical  region  the  frequency  has  just  reached  1/2 /rev, 
and  the  root  locus  is  at  the  point  A.  For  still  larger  V 
the  frequency  remains  fixed  while  the  real  part  of  one 
root  decreases  and  that  of  the  other  Increases.  When 
y reaches  the  value  for  which  the  hover  root  has  a fre- 
quency of  1 /2/rev,  the  locus  is  at  the  center  of  the 
critical  region;  there  the  roots  have  their  maximum 
stability  change  (which  Is  0(p) ) so  the  locus  is  at  the 
point  B.  As  > increases  more,  the  locus  moves 
toward  the  other  boundary  of  the  critical  region.  The 
locus  reaches  that  boundary  at  the  point  C , and  for 
still  larger  y the  frequency  is  no  longer  fixed  at  1/2/ 
rev;  rather  the  real  part  of  the  root  is  the  same  as  the 
hover  value,  while  there  is  an  0(p>  change  in  the 
frequency  which  decreases  in  size  as  y increases. 

For  still  larger  y the  locus  is  again  identical  (to  order 
p)  to  the  hover  locus. 

Figure  2 shows  typical  root  loci  for  fixed  T »nd 
vaiying  p . The  circle  the  locus  starts  from  is  t'.e  y 
locus  for  hover  (p  = 0)  and  the  appropriate  Kp  The 
y for  each  locus  may  be  found  from  Be  A at  p = «.  , 
since  for  the  hover  root  Re  A = - y/16  . As  p in- 
creases from  zero,  for  the  roots  near  ImA  =1/2 
there  Is  an  0<p)  change  In  the  frequency  pulling  the 
root  toward  ImA  = l/2  ; while  the  damping  remains 
fixed  at  the  hover  value  The  locus  reaches  ImA  = l/2 
for  a corner  p which  corresponds  to  the  boundary  of 
the  critical  region.  For  larger  p , frequency  remains 
lixed  at  1/2 /rev  while  one  branch  of  the  locus  moves  to 
the  left  (decreased  stability)  and  the  other  to  the  right 
(Increased  stability).  This  behavior  of  the  p loci  is 
characteristic  of  roots  of  a system  with  periodic 
coefficients. 

2 

Order  p Results 

o 

To  order  p the  roots  for  most  * and  yQ  are 


* “ “ T6  + 1 { 1,2  * (1  + 


-II 


II) 

V 18/ 


2 8 


2 y 

’ -8KP+‘,K, 


9 /y\2  2 y i. 

\i l)  ~v  -tVi1 


(8) 


and  its  conjugate.  Thus  away  from  ImA  = l/2  or  1 the 
roots  are  just  the  hover  roots  with  an  0<p2)  change  in 
the  frequency.  There  are  two  effects  ot  p ; the  first 
corrects  the  term  y/8  Kp  in  the  hover  frequency  to 
properly  account  for  the  average  of  Kp  M s . The 
second  effect,  that  in  the  last  term  of  the  frequency,  is 
ntlrely  due  to  the  periodic  aerodynamic  woeti/cients; 
this  is  the  first  effect  of  the  periodic  coefficients  seen 
in  the  analysis,  except  for  the  critical  region  near 
1mA  * l/2.  Typical  root  loci  for  varying  p are 
shown  in  Fig.  2.  These  are  the  loci  that  are  not  near 
ImA  = l/2  or  1 ; the  frequency  change  is  small  even 
at  p = 0.5  . Eq.  8 may  also  be  used  for  the  branches 
of  the  root  loci  on  the  real  axis  when  the  quantity  under 
tb  jquare  root  sign  is  negative  (i.e. , for  y large 
e>;  igh).  There  are  two  real  roots  then,  the  (+1)  in  the 
f’  equency  becoming  (*1) . A point  on  the  locus  of 
special  interest  Is  where  one  branch  of  the  Iocub  on  the 
real  axis  crosses  into  the  RHP.  i.e. , becomes  un- 
stable (as  in  Fig.  1).  The  criterion  for  this  divergence 
boundary  is  A = o , or  to  order  p2  ; 
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The  effect  of  p on  the  RHS  (due  to  the  periodic  co- 
efficients) dominates  mat  on  the  LHS  (due  to  the  average 
of  KpMd)  for  all  values  of  y and  v . Thus  the  critical 
value  of  negative  Kp  , beyond  which  the  locus  lies  in 
the  RHP  , is  actually  increased  by  increasing  p . The 
criterion  from  the  hover  case  is  conservative  then; 
this  is  the  opposite  of  the  conclusion  that  would  have 
been  reached  from  a consideration  of  the  averaged  co- 
efficients only. 

2 

To  order  p there  is  a critical  region  when 
ImAg  = 1 , I.e. , when 


The  behavior  of  the  loci  near  the  critical  region  Is 
similar  to  that  near  ImAg  = l/2  , except  that  now 
changes  are  0(p2)  , not  0(p)  . The  boundary  of  the 
critical  region  is  described  by  a relation  between  v2 
and  T2  with  •’i  “ = 0 ; the  critical  region  is  thus 
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a narrow  band,  of  width  0<p2) , about  Ira  A 1-.  In- 
side the  critical  region  the  frequency  is  fixed  at  ImA  « l 
with  an  O (p2)  change  in  the  damping.  The  maximum 
reduction  ’•*  enhancement)  of  the  damping  Is 


o 

which  is  an  0(p  ) small  reduction  in  the  large  hover 
damping.  Figs.  1 and  2 show  typical  root  loci  near 
ImA  = 1 for  varying  Y and  p . 

Flap  Hate  Feedback 

The  use  of  flap  rate  feedback,  KR  / 0 . does  not 
change  the  behavior  of  the  solution  quantitatively.  The 
hover  root  becomes 
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and  there  are  critical  regions  about  ImAg  =1/2  and  1 
again.  The  critical  region  boundaries  and  stability  de- 
gradation depend  on  KR  now.  It  is  necessary  that 
Kr  > -1  for  the  hover  root  to  be  stable,  but  Kr  > 0 
will  be  the  usual  case  anyway. 

Y-p  Plane 


(one  more  and  one  less  stable  than  the  hover  root) . This 
behavior  is  Just  that  seen  already  in  the  p loci  (Fig.  2). 
Figs.  3,  4,  and  5 may  be  compared  with  similar  ones  In 
Ref.  2,  which  were  constructed  from  numerical  calcula- 
tions; on  the  basis  of  this  comparison,  the  O (p2)  ana- 
lytic results  are  quite  accurate  up  to  p * 0.5  or  so. 
There  is  some  discrepancy  between  the  results  for  the 
ImA  = 1 region  however,  particularly  with  Kp  = -0.1  , 
although  the  change  in  scale  (Ref.  2 shows  results  out  to 
P = 2.5)  exaggerates  the  difference.  For  v exactly  1 
the  analytic  results  indicate  no  critical  region  If  Kp  < 0; 
but  only  a slightly  larger  (for  example  m = 1.01)  is 
necessary  to  get  a sizable  critical  region  with  Kp  = -0. 1 
(see  Fig.  5).  The  analytic  results  show  the  v = 1 case 
is  a very  sensitive  one  for  small  Y and  small  Kp  , 
and  it  is  unlikely  that  a numerical  calculation  could  treat 
the  case  accurately.  Of  course  an  actual  rotor  will  al- 
ways have  v at  least  slightly  greater  than  1 , so  the 
numerical  calculations  are  probably  reliable  then; 
furthermore,  the  discrepancy  may  also  be  an  indication 
that  for  very  small  Y the  analytic  results  (to  order  p2) 
are  not  valid  out  to  as  large  a p as  they  are  for  more 
reasonable  Y . In  any  case  this  discussion  illustrates 
the  kinds  of  problems  thv.  may  be  hidden  in  a purely 
numerical  solution;  they  can  only  be  found  and  studied 
by  analytic  procedures  (which  at  least  tell  where  to  look 
for  problems). 

The  Small  Y Case 

The  small  parameter  in  this  case  is  the  Lock  num- 
ber Y . Setting  Y = 0 removes  all  the  aerodynamic 
terms,  and  the  roots  are  A = tip  for  all  p . The 
method  of  multiple  time  scales  is  used  to  obtain  solu- 
tions uniformly  valid  for  all  # . The  time  scales  con- 
sidered are  of  order  1,  of  order  yl  , etc.  P ij 
expanded  as  a series  In  Y (P  = Pg  + YPj  + . . .)  and 
sols  » (v  - vQ  + Ypj  ♦ . . .).  The  advance  ratio 
p is  not  small  now,  and  the  reverse  flow  r"?‘ — .nusi  Le 
considered. 


The  results  of  the  small  p analysis  may  be  used  to 
plot  lines  jf  constant  Re  A and  ImA  on  the  Y-p  plane. 
Typical  results  are  shown  in  Figs.  3,4,  and  5 for 
v - l.o  and  Kp  = 0,  0.1  , and -0.1,  respectively. 

The  critical  regions  appear  in  the  Y-p  plane  as  regions 
in  which  ImA  is  constant  (at  l/2/rev  or  1/rev);  they 
are  indicated  in  the  figures  by  the  circled  values  of 
ImA  (the  region  where  Im>  = 0 is  where  there  are  two 
real  roots,  not  a critical  region).  These  figures  are 
interpreted  as  follows.  A horizontal  line  is  a line  of 
constant  Y , and  so  as  p varies  It  gives  the  corre- 
sponding values  of  Re  A and  ImA  just  as  a p root  locus 
does.  Similarly  a vertical  line  is  a constant  p line, 
and  so  gives  A as  a function  of  Y just  as  a Y root 
locus  does.  For  example,  consider  a horizontal  line  in 
Fig.  3 (Kp  = 0)  with  Y = 8.  i.e. , Y/i6  = 0.5.  As 
p increases,  the  line  remains  parallel  to  the  Re  A = 
constant  lines  so  ReA  remains  fixed  at  the  hover  value. 
The  ImA  = 1/2  region  comes  closer  to  the  horizontal 
line  as  p increases,  which  means  that  ImA  moves 
toward  1/2 /rev.  Eventually  the  constant  Y line  crosses 
into  the  ImA  = 1/2  region;  then  ImA  is  fixed  at  1/2/ 
rev  while  for  each  point  In  the  region  there  are  two 
values  of  ReA , giving  the  damping  for  the  two  branches 


Order  Y Results 

If  v 4 n/2  for  any  integer  n (i.e.  ImAQ  not 
equal  to  a multiple  of  1 /2/rev)  the  root  to  O(y)  If 
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where  Mp  , Mp  and  Mg  are  the  averages  of  the  aero- 
dynamic coefficients  around  the  azimuth.  For  the  co- 
efficients considered  here,  Mp  = 0 for  all  p . Flap 
rate  feedback  only  effects  (to  order  Y)  the  damping;  the 
real  part  of  may  be  written 


ReA 


■(1) 


Mi 


0 . 0 
-Mp/Mg> 


(14) 


The  ratio  -M^/M®  determines  the  relative  effect  of 
Kr  ; this  ratio  is  shown  In  Fig.  6 as  a function  of  p . 
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This  parameter  la  a positive  number,  wu'ch  varlea  'lttle 
with  p (from  1 a*  p = 0 to  1 /2  at  p = • , with  moat 
of  the  change  below  p » 1) . The  negative  of  this  ratio 
gives  a critical  KR,  since  it  is  necessary  that 

KR  > Kr  * M?  M in  order  that  the  system  be 
crit  ^ 

stable  for  small  Y (to  O(y)  and  with  t>  4 n/2) . 

For  K^  = 0 the  root  is 
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+ it' 
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The  lerodynamlc  coefficients  -sM  and  8 Mg  are 
alwa  m positive;  they  have  the  value  1 for  p = 0 and  are 
asymptotic  to  (8/3r)p  and  (16/3*)p , respectively, 
for  large  p ; these  coefficients  are  shown  in  Fig.  7. 

The  root  loci  for  varying  Y and  varying  p are 
shown  in  Fig.  8 for  Kp  = 1 , KR  - 0 , and  v = l ; 
the  locus  for  Kp  = - 1 is  obtained  by  reflecting  this 
locus  about  the  ImX  ■ v line.  The  locus  for  Kp  = P 
is  difficult  to  plot  since  ImX  = v for  all  Y and  p (to 
order  Y) ; it  may  be  visualized  by  projecting  the 
Kp  - 1 loci  onto  theline  ImX  - v.  These  loci  should 
be  compared  with  the  small  Y portions  of  the  curves  in 
Fig.  1,  which  are  for  small  p . In  Fig.  8.  the  Y locus 
for  a given  p starts  out  at  X •*  iv  always,  and  is  a 
straight  line  with  slope 


0ImX 
0ReX 
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which  varies  from  -Kp/v  to  -2Kp/v  for  p from 
0 to  «>  (see  Fig.  6 for  -m2  /M°.  ) . The  step  sire  on 
the  Y locus,  for  a unit  change  in  Y/lfi  , is 


tor  p from  0 to  1 to  *>,  respectively.  The  p locus 
for  a given  Y starts  out  virtually  from  the  p = 0 line 
and  is  asymptotic  to  the  p ■ » line,  with  the  step  size 
on  the  locus  for  a unit  chaise  in  p increasing  as  p in- 
creases. For  reasonable  p the  locus  does  not  vary 
much  from  the  small  p results.  An  O(Y)  analysis 
can  only  obtain  the  slope  of  the  Y locus  at  Y = 0 , so 
the  locus  is  a straight  line,  as  found  above;  to  find  the 
curvature  effect  it  Is  necessary  to  go  to  order  Y2  . The 
significance  of  the  curvature  (order  Y2)  may  be  judged 
from  a comparison  of  the  Y loci  of  Fig.  8 (all  p . small 
Y)  and  Fig.  1 (small  p , all  Y) ; on  this  basis  the 
small  y results  should  be  limited  to  Y/i6  less  than 
0.2.  On  the  basis  of  neglect  of  the  curvature  effects 
alone  the  results  might  be  accepted  to  higher  Y , but 
the  O(Y)  results  will  also  be  limited  by  the  effects  of 
the  critical  region,  which  will  be  examined  next. 

The  discussion  of  the  O(Y)  analysis  has  so  far 
only  been  concerned  with  the  basic  roots,  meaning  the 
roots  away  from  the  influence  of  a critical  region.  In 
the  small  Y case,  the  criterion  for  being  away  from  a 
critical  region  is  that  4 (n/2)  for  any  integer  n ; 
this  may  be  written 


v 


i.e. , the  rotating  natural  frequency  may  not  be  a distance 
of  order  Y/16  from  n/2/rev.  Since  v is  almost  al- 
ways just  above  1/rev  this  criterion  is  seldom  fulfilled, 
and  the  critical  regions  may  be  expected  to  dominate  the 
root  loci  behavior  for  small  Y . Furthermore,  if  Kp 
is  large  enough  positive  or  negative,  the  basic  locus  will 
also  cross  ImX  = 3/2  or  l/2  for  Y/16  still  small 
(see  Fig.  8),  so  these  critical  regions  may  affect  the 
loci  even  if  the  ImX  = 1 region  does  not. 

If  = n/2  for  some  integer  n , there  arise 
critical  regions,  with  behavior  of  the  root  loci  similar 
to  that  encountered  already  in  the  small  p case.  The 
case  Kr  = 0 will  be  considered  first.  The  critical 
region  boundary  is  given  by 


which  varies  from 
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where 
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and  similarly  for  and  Mg  (i.e.,  the  nth  harmonic 
in  a complex  Fourier  series  representaMon  of  the 
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aerodynamic  coefficients).  This  equation  describes  a 
band  of  width  O(Y)  about  ImA  « n/2 . The  maximum 
stability  chance  occurs  at  ►'j  * - (Kp/2vo)  Mj  , that  la 
at  the  center  of  the  critical  region,  where  the  basic  root 
(Eq.  15)  would  cross  1mA  ■ t'o  m n/2  . There  the  root 
is 
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Outside  the  critical  region  there  is  an  order  Y change 
in  the  frequency  of  the  basic  root  with  no  change  In  the 
damping.  Inside  the  critical  region  the  frequency  is 
fixed  at  n/2 /rev  while  there  is  an  older  Y change  In 
the  damping.  The  damping  of  the  basic  root  is  itself 
O(Y)  however,  so  In  contrast  to  the  small  it  case,  the 
critical  region  can  here  lead  to  actual  Instability,  not 
just  stability  degradation. 


from  the  basic  root  damping,  (V/2)Mp.  i a always  nega- 
tive (Fig.  7);  as  for  the  small  it  case  it  is  O(y)  , but 
here  that  means  the  basic  damping  is  amall.  In  fact  it 
la  the  same  order  as  the  critical  region  contribution,  so 
the  destabilised  root  may  be  actually  unstable,  rather 
than  just  a small  perturbation  from  the  basic  damping 
as  for  the  small  p case.  The  behavior  of  the  locus 
depends  on  the  relative  effects  of  M?  and  the  nth  har- 
monica of  the  aerodynamic  coefficients  under  the  square 
root  in  Eq.  (17).  For  most  cases  the  critical  region 
effect  dominates,  so  that  as  p increases  it  eventually 
reaches  a critical  value,  at  which  point  (for  the  case  of 
maximum  stability  change)  one  root  crosses  into  the 
RHP,  i.e.,  becomes  unstable.  From  Eq.  17,  it  follows 
that  increasing  Kp  (Kp  either  positive  or  negative) 
always  increases  the  effect  of  the  criticai  region  which 
means  decreasing  the  critical  p for  which  the  root 
becomes  unstable.  Thus  the  critical  p is  a function  of 
I Kp  | , for  each  of  the  critical  regions;  this  function 
may  be  found  from  Eq.  (17)  by  setting  ReA  *=  0 (the 
requirement  for  crossing  the  1mA  axis).  Since  the 
aerodynamic  coefficients  are  rather  complex  functions 
of  p , it  is  more  convenient  to  find  the  critical  | Kpl 
as  a function  of  p : 


The  root  loci  for  small  Y show  the  behavior  char- 
acteristic of  periodic  systems,  and  familiar  from  the 
discussion  of  the  loci  for  the  small  p case.  For  some 
p or  Y the  locus  will  cross  into  the  critical  region. 

The  maximum  stability  change  inside  the  critical  region 
occurs  at  the  center,  which  Is  where  the  basic  root 
locus  (e.g. , as  shewn  in  Fig.  8)  would  have  crossed 
1mA  - n/2  . Whether  the  center  is  ever  reached  de- 
pends on  the  parameters  involved  (Kp  , v , Y , p ) . For 
the  O(Y)  analysis,  the  locus  has  a slope  proportional 
'jo  Kp  , so  for  | Kp  | not  too  small  the  locus  will 
probably  cross  some  ImA  ■ n/2  line  with  V/16  still 
small  (e.g. , in  Fig.  8 where  Kp  - 1 , the  loci  cross 
the  ImA  - 3/2  line).  If  Kp  » 0 however,  there  is  no 
O(Y)  change  in  the  frequency,  so  ImA  - v for  all  p 
and  Y . Then  the  locus  will  never  reach  the  center  of 
the  critical  region  (cross  ImA  ■ n/2)  unless  v » n/2  , 
in  which  ctse  the  locus  is  always  at  the  center.  More 
generally,  if  |Kp|  is  very  small,  the  locus  will  not 
cross  ImA  » n/2  until  Y/l6  is  so  large  as  to  be  out- 
side the  range  of  the  order  Y analysis.  An  0(Y  ) 
analysis  would  change  significantly  the  conclusions  about 
whether  the  center  of  the  critical  region  is  reached  under 
certain  conditions,  particularly  for  small  | Kp  | when 
the  0(Y2)  change  in  the  frequency  is  more  important 
than  the  O(Y)  change.  For  example,  Fig.  1 Indicates 
that  with  v slightly  greater  than  1 the  Y locus  (for 
small  it  ) would  never  cross  ImA  = 3/2  for  Kp  * 1 
while  it  would  always  cross  ImA  = 1 for  Kp  ■ 0 , just 
the  opposite  of  the  conclusions  Indicated  by  the  O(Y) 
results  (Fig.  8).  In  any  case,  since  the  maximum 
stability  change  occurs  at  the  center  of  the  critical 
region,  it  is  useful  to  examine  it  as  a worst  possible 
case,  which  may  perhaps  be  approached  but  never 
reached  for  certain  values  of  v and  Kp  . 

Considering  the  p root  locus  again , the  maximum 
change  in  the  damping  from  the  value  of  the  basic  root 
is  given  by  the  real  part  of  Eq.  17.  The  contribution 


(18) 


This  may  be  regarded  as  a maximum  | K p I for  a given 
p ; for  larger  | Kp  I the  locus  is  in  the  RHP  at  that  p . 
These  boundaries  of  |Kp|  vs.  p are  shown  in  Fig.  9 
for  ImA  * 1/2  , 1 , and  3/2  . With  the  exception  of 
the  roots  near  ImA  • 1/2  with  p above  0.5  or  so, 

Fig.  9 shows  the  criterion  on  |Kpj  is  not  very 
stringent.  Fig.  9 shows  also  that  near  ImA  * l the 
roots  are  always  stable,  regardless  of  p , if 
I Kp  | < Ji  ; the  loci  may  be  expected  to  be  near 
ImA  =1  for  zero  or  small  Kp.  In  terms  of  the  p 
locus , this  means  that  as  p increases  the  locus  does 
not  cross  into  the  RHP  . Just  after  the  locus  crobttes 
the  critical  region  boundary,  the  effect  of  the  criticai 
region  is  seen  and  one  branch  moves  to  the  right  and 
the  other  to  the  left  (as  do  the  loci  in  Fig.  2).  As  p 
increases  further  however,  the  lamping  of  the  basic 
root  (which  is  always  stable,  and  increases  with  p ) 
eventually  dominates  the'effect  of  the  critical  region, 
and  the  root  which  wa?  becoming  less  stable  turns 
around  before  reaching  t’e  RHP  . So  for  still  larger  p 
both  branches  of  the  locus  ."ill  be  moving  to  the  left, 
i.e.  becoming  more  stable  as  p increases.  For  the 
root  loci  near  ImA  = 1/2  or  3/2  , the  effect  of  the 
critical  region- remains  dominant,  and  so  one  root 
eventually  crosses  into  the  RHP  as  p is  increased. 

The  critical  p is  considerably  lower  for  ImA  = 1/2 
than  for  ImA  - 3/2  . This  points  out  an  undesirable 
feature  of  negative  pitch  flap  coupling,  Kp  < 0 ; not  so 
much  that  it  reduces  the  critical  p , but  rather  that  it 
moves  the  basic  root  nearer  to  ImA  = 1 /2  . 
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The  om  of  flap  rate  feedback,  KR  i 0 , results  In 
(x>  qualitative  changes  In  the  behavior  of  the  loci.  KR 
la  however  a useful  design  parameter;  It  may  be  used 
for  example  to  raise  the  critical  p or  I Kp  | . 

Evaluation  of  the  Order  Y Results 


Numerical  calculations  were  made  of  the  p root 
loci  for  moderate  and  small  values  of  Y . CX,  the  basis 
of  a comparison  of  the  numerical  and  analytic  results . 

It  Is  concluded  that  the  small  Y analysis  to  order  Y is 
useful  only  for  truly  small  Y , t.g. . Y = 2 or  3 
(Y/16  = 0.2  or  so).  Problems  are  encountered  with 
both  the  basic  roots  and  the  effects  of  the  critical  region . 
The  basic  root  to  order  Y neglects  the  curvature  of  the 
Y locus,  which  is  especially  Important  for  zero  or 
small  Kp  , since  then  the  change  of  ImX  for  small  y 
is  due  more  to  0(Y2)  terms  than  to  the  O(Y)  term. 
The  damping  of  basic  root  Is  O(Y)  always,  no  matter 
what  order  the  analysis  is  carried  too;  for  example  the 
0(p2)  results  give  ReX  «-  Y/ie  for  all  Y.  However 
while  basic  damping  is  O(Y).  the  contribution  to  the 
damping  due  to  the  critical  region  will  have  terms  that 
are  0(Y2).  Thus  for  large  enough  Y the  conclusions 
in  the  dl*eussion  above  of  the  effects  of  the  critical 
region  on  the  p root  loci  will  not  be  valid,  since  they 
depend  on  the  basic  and  critical  region  damping  being  of 
the  same  order  in  Y . In  particular,  the  behavior  of  the 
locus  in  which  the  root  being  destabilized  by  the  critical 
region  turns  around  and  becomes  more  stable  due  to  the 
eventual  dominance  of  the  basic  damping  is  not  possible 
except  for  very  small  Y.  for  which  0<y2)  effects  are 
ii.  fact  negligible.  Indeed,  it  was  found  in  the  numerical 
calculations  that  with  Y = 6 (Y/16  * 0.375,  which  is 
not  very  small),  v near  1,  and  Kp  zero  or  small  so 
the  root  is  near  1/rev,  that  the  p locus  does  not  turn 
around  but  rather  eventually  crosses  into  the  RHP.  The 
stability  boundaries  given  in  Fig.  9 are  only  valid  then 
for  truly  small  values  of  Y/16  . 


Order  Y2  Results 


To  order  Y‘  the  basic  root,  valid  for 
is  (with  Kr  *=  0) 


t n/2  . 
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As  reported  above,  there  is  no  0(Y  ) changeJn  ReX  . 
and  the  0(Y2)  change  in  the  frequency  is  dominant  for 
small  Kp  . indeed  for  Kp  - 0 the  only  change  in 
ImX  is  0(y2).  The  0(Y2)  results  would  significantly 
alter  the  plots  of  the  basic  root  loci  shown  in  Fig.  8. 


The  Large  Y Case 

The  small  parameter  In  this  case  is  the  Inverse  of 
the  Lock  number.  T'1  . For  Y very  large,  the  aero- 
dynamics dominate  the  system.  For  Y «=  « the  inertia 
and  centrifugal  spring  terms  are  negligible,  leaving  a 


first-order  system  that  does  not  depend  on  Y . This 
system  gives  one  real  root,  which  Is  independent  of 
and  valid  for  Y * •• . The  reduction  of  order  of  the 
equation  at  Y = • Isa  characteristic  of  a boundary- 
layer  type  of  problem.  This  characteristic  means  that 
one  of  the  roots  goes  to  - * as  Y goes  to  •» . The  two 
roots  may  be  found  by  use  of  the  following  substitution. 


f 1 

P - exp  f y p_j  ♦ p„  + y p_j  ♦ . . . d* 

This  substitution  gives  two  solutions,  from  which  the  two 
roots  are  obtained  as  (for  K^  » 0) : 
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Thus  there  are  two  real  roots,  one  (X2)  approach- 
ing - «o  as  Y increases  to  » (Mp  is  negative) . and  the 
other  (>  ) approaching  a constant.  This  behavior  of 
the  Y • oot  loci  is  expected  from  the  small  p results  ; 
Fig.  1 shows  that  for  large  enough  Y the  locus  is  on  the 
real  axis,  i.e.  there  are  two  real  roots,  one  approaching 
- « and  the  other  -Kp  for  Y — » . To  lowest  order 
Xj  does  not  depend  on  Y . because  it  represents  the 
balance  of  the  aerodynamic  damping  and  the  aerodynamic 
spring  only.  The  value  of  Xj/(-Kp)  for  varying  p . and 
Y = ««  . is  shown  in  Fig.  10;  the  movement  show",  takes 
place  entirely  on  the  real  axis  in  the  X plane.  As  for 
the  small  g case  (Fig.  1)  the  root  on  the  real  axis  is  in 
the  LHP  if  Kp  > 0 and  in  the  RHP  - unstable  - if 
Kp  < 0 . The  value  of  Xj/(-Kp)  varies  from  1 to  7 /8 
for  p = 0 to  • , with  most  of  the  change  between 
p = 0.5  and  g = 1 ; thus  there  is  little  variation  of  the 
root  with  g (to  order  1 in  Y).  The  size  of  the  0(Y_1) 
term  in  X j is  indicated  by  the  result  for  p =0  which 
is  easily  obtained  (since  the  aerodynamic  coefficients  are 
constant  then)  as 
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To  lowest  order  A2  Is  - (Y/16)  2^-8  Ma]  . The 

0 ^ 

aerodynamic  coefficient  -8Mfl  is  given  in  Fig.  7.  For 

4 — 1 it  has  the  value  -8M^  « 1 + p'/8;  for  large 
4 it  is  asymptotic  to  (8/3*)  4.  This  root  becomes  in- 
creasingly negative  as  Y Increases,  and  also  as  4 
increases.  The  order  1 term  in  X9  is  the  negative  of 
the  lowest  order  term  in  Xj ; thus“the  behavior  of  this 
term  is  also  given  by  Fig.  10. 

The  solution  for  the  roots  is  not  valid  if  Mjj  = 0 for 
any  ^ . When  Mp  is  small,  the  assumptions  made 
about  the  order  of  terms  in  deriving  the  solutions  are 
violated;  this  maybe  inferred  from  the  expressions  ob- 
tained for  the  roots  by  the  continual  appearance  of  Mp 
in  the  denominator  of  the  integrands  (Eqs.  (20  and  (21) . 
There  is  a narrow  transition  region,  of  width  0(Y-^'^)  . 
about  any  points  where  Mp  = 0 , i.e.  , where  the 
damping  goes  through  zero.  As  it  happens,  however, 

Mp(  (1 . 4)  is  a negative  quantity  which  never  reaches 
zero.  Thus  for  KR  = 0 there  arc  no  transition  regions 
or  boundary  layers,  and  the  solutions  obtained  from  the 
substitution  for  (3  are  uniformly  valid  over  the  entire 
azimuth. 

Flap  Rate  Feedback 

When  Kn  = 0 , there  are  no  boundary  layers  or 
transition  regions  because  Mp  < 0 always.  With  flap 
rate  feedback,  KR  / 0 , the  same  expressions  for  the 
main  solutions  for  (1  are  obtained  as  for  = 0 , ex- 
cept that  Mjj  is  replaced  by  Mp  - KR  . The  aero- 
dynamic coefficient  - (Mp  - KR  M^)  can  become  nega- 
tive over  regions  of  the  disk  for  certain  combinations  of 
4 and  Kr;  that  is,  there  muj  be  negative  damping 
over  part  of  the  azimuth  range.  When  such  regions  of 
negative  damping  exist  it  means  there  must  be  transition 
regions  about  the  points  where  the  damping  goes  through 
zero.  Since  Eqs.  (20)  and  (21)  for  Xj  and  X2  were 
derived  on  the  basis  of  no  transition  regions,  those  ex- 
pressions for  the  roots  (with  Mp  replaced  b/  Mp  - 
KRMfl)  are  not  vaiid  when  there  is  negative  damping 
over  any  part  of  the  disk. 

The  criterion  for  the  existence  of  transition  regions 
is  that  there  be  negative  damping  on  some  portion  of  the 
disk,  i.e.,  - (Mp  - 0.  Mp  is  always 

negative;  is  usually  positive,  but  may  be  negative 
on  the  retreating  side  for  large  enough  4 . If  Kpj  is 
too  large  positive,  the  negative  values  of  Mg  on  the 
retreating  side  eventually  dominate  Mp  as  4 is  In- 
creased, so  there  will  be  negative  damping  on  the  re- 
treating side;  if  is  too  large  negative,  KrM# 
eventually  dominates  Mp  on  the  advancing  side  and 
there  wiii  be  negative  damping  there  if  4 is  large 
enough.  Quantitative  values  of  maximum  and  minimum 

as  a function  of  4 are  given  in  Fig.  11.  For  the 
cases  with  negative  damping  there  will  be  transition 
regions  (of  width  0(Y“^^))  near  where  Mp-K^M^ 

0 which  greatly  complicates  the  analysis.  For  these 
cases  it  is  also  expected  that  there  wiil  be  other  prob- 
lems, including  numerical  computation  problems, 
physical  control  problems,  and  large  flapping  ampli- 
tudes. Thus  while  a region  of  negative  damping  doe3 
not  necessarily  mean  there  is  a Happing  instability,  it 


does  mean  that  there  are  many  problems  - ana- 
lytical, computational,  and  physical  - so  requiring 
-(Mp  - Kp M ^ ) > 0 is  a reasonable  design  criterion. 
This  criterion  provides  a maximum  and  minimum 
for  a given  4 . The  limits  of  from  this  rue  are 
much  easier  to  obtain  than  actual  stability  boundaries; 
and  Fig.  11  shows  that  although  conservative,  it  is  not 
a serious  restriction  for  4 less  than  1 or  2.  For  large 
4 It  is  a serious  limitation,  indicating  that  M#  (flap 
moment  due  to  blade  pitch)  is  not  very  good  for  flapping 
rate  feedback  then.  Although  the  derivation  of  this  ruie 
has  been  based  on  the  iarge  Y case,  the  criterion  of  no 
negative  damping  has  nothing  to  do  with  Y , and  so 
should  be  a reasonable  criterion  for  ail  Y.  Indeed  the 
criterion  > -1  for  4 0 is  the  same  as  from  the 

smaii  4 case,  where  it  is  a true  stability  criterion, 
and  valid  for  all  Y . 

The  Large  4 Case 

The  small  parameter  for  this  case  is  the  inverse  of 
the  advance  ratio,  4“*  . For  4 very  large,  the  aero- 
dynamics again  dominate  the  system.  However,  for 
Kft  = 0 the  damping  is  Ofcj)  while  the  aerodynamic 
spring  is  O (g2)  . with  the  result  that  at  4 — • » the 
aerodynamic  spring  must  be  balanced  by  the  inertiai 
forces  in  order  to  obtain  an  equation  with  the  proper 
order  of  terms.  Thus  the  equation  to  lowest  order  in 
4 takes  the  form 


p ♦ 4C(*)P  + 42K(i l>)  P = 0 (22) 

The  solution  of  this  equation  is  either  a rapid  sinusoidal 
osciitation  with  frequency  of  0(g  ) , or  a sum  of  ex- 
ponentials with  time  constants  of  0(g"l)  , depending 
on  whether  the  aerodynamic  spring  is  negative  or  posi- 
tive (the  criterion  is  a bit  more  complicated  reaiiy , but 
that  statement  wiii  dc  for  the  present  discussion).  The 
aerodynamic  spring  changes  sign  in  the  middie  of  the 
advancing  side  and  again  in  the  middle  of  the  retreating 
side,  and  at  each  point  there  is  a transition  region  (of 
width  O (g-^/2)  ) across  which  the  solutions  must  be 
matched.  There  arc  also  transition  regions  (of  width 
0(4"2/3)  j between  the  advancing  and  retreating  sides 
of  the  disk  (near  0 « 0 and  180°)  through  which  the 
solution  must  be  matched;  near  these  points  the  aero- 
dynamic spring  is  very  smaii,  although  it  does  not 
change  sign.  In  general,  near  any  point  where  the  aero- 
dynamic spring  is  very  small,  the  assumption  that  aii 
terms  in  Eq.  22  are  of  the  same  order  is  violated,  and 
there  must  be  a narrow  transition  region  about  such  a 
point.  In  this  problem  there  are  four  such  points, 
dividing  the  disk  into  four  regions,  each  with  its  own 
main  solution.  These  main  solutions  may  be  obtained 
using  the  substitution 


This  substitution  gives  the  main  solution  (to  order  p ) 
as  follows. 


10 


„ , |(|-T  T5*' 


c 


(23) 


where 


f = j I sin  0 | (-^  | sin  ^ | - cos  $ - Kpsintf> ) 

g = 1 ( {cos*  - Kp  sin  ip  - ■—  sin  ip ) 

( * 1 on  the  advancing  side 

r = 

(-1  on  the  retreating  side 


and  Cj  and  C2  are  constants  (If  f < 0 then  Cj  is 
complex  and  C2  = Cj ) . There  is  a main  solution  of 
this  form  for  each  of  the  four  outer  regions  of  the  disk; 
the  solutions  are  matched  through  the  four  transition 
regions  (giving  the  two  constants  on  one  side  in  terms  of 
the  two  on  the  other  side)  by  the  standard  techniques  of 
perturbation  theory.3,5  By  this  procedure  the  transient 
solution  for  given  initial  conditions  may  be  constructed 
around  the  entire  disk.  Floquet  theory  then  gives  a 
quadratic  equation  for  the  two  roots  as  follows. 
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where  b is  a function  of  p , V , and  Kp  . The  solu- 
tion for  the  roots  is  the  i 


b + ni  for  b > 1 

b + ni  for  -1  < b < 1 

| b I + — + ni  for  b < -1 


where  n is  some  integer.  This  result  shows  the 
typical  behavior  of  the  roots  of  periodic  svstems.  For 
Ibi  < 1 the  damping  is  fixed  at  -p  (T/Gw)  with  a 
change  due  to  b in  the  frequency;  for  b > 1 the  fre- 
quency is  fixed  at  n/rev  with  a positive  and  negative 


change  due  to  b in  the  damping;  for  b < -1  the 
frequency  is  fixed  at  n + 1 /2/rev  with  a positive  and 
negative  change  in  the  damping.  The  critical  region 
boundaries  are  given  by  b = 1 and  b = -1 . 

The  general  character  of  the  critical  regions  and 
instability  boundaries  in  the  V - p plane,  as  obtained 
from  the  solution  of  Eq.  (24),  is  sketched  in  Fig  12. 
Because  p is  large , it  happens  that  | b | is  much 
greater  than  1 almc  n ilways,  so  the  critical  regions 
dominate  the  behav.„.  of  the  roots.  The  sign  of  b 
changes  regularly  however;  b must  of  course  go  through 
zero  then,  but  it  does  so  very  quickly,  so  there  is  only 
a very  narrow  band  between  the  ImX  n/rev  and  the 
ImX  = I.  + l/2/rev  regions  in  which  I b [ < 1.  When 

Ibl  < 1,  the  real  part  of  X is  -p  (T/Cir)  , i.c.,  the 
root  is  stable  for  all  p and  V ; thus  there  must  always 
be  a band  of  stability  surrounding  the  transition  from 
n/rev  to  n + l/2/rev.  These  characteristics  are 
illustrated  in  Fig.  12.  The  locus  between  the  critical 
regions  has  a rather  fine  structure  which  would  be 
difficult  to  obtain  numerically.  A root  locus  for  varying 
y or  p (a  vertical  or  horizontal  section  in  Fig.  12)  in 
the  vicinity  of  a critical  region  boundary  would  in  quick 
succession  move  from  the  RHP  (unstable)  to  the  LHP 
(stable)  with  frequency  fixed  at  n/rev,  rapidly  move 
from  ImX  = n/rev  to  ImX  = n + 1/2/rev  in  the 
RHP  with  damping  given  by  -p  tV/Gir)  (which  would 
be  nearly  constant  because  the  critical  region  boundaries 
are  so  close),  and  then  move  from  the  LHP  into  the  RHP 
with  frequency  fixed  at  n + 1/2 /rev. 

Fig.  12  shows  that  for  a given  p the  system  Is 
stable  for  a large  enough  V . Positive  Kp  is  stabiliz- 
ing, tending  to  decrease  the  size  of  the  instability  regions; 
negative  Kp  is  destabilizing  in  this  sense.  The  rotating 
natural  frequency  of  the  flap  motion,  v , docs  not  enter 
the  high  p case  to  order  pg  (the  acrodynam.c  soring 
dominates  the  centrifugal  spring  until  order  pj) ; this 
is  consistent  with  the  fact  the  critical  regions  dominate 
the  high  p behavior,  so  the  frequency  is  fixed  at  a 
multiple  of  1/2/rev. 


A comparison  of  these  analytical  rcsuMs  with  the 
results  of  numerical  calculations  indicates  that  the  high 
p solution  is  good  down  to  p =2.5  or  so.  Thus 
numerical  calculations  are  required  to  join  the  loci 
from  p = 0.5  to  2.5  say  (for  V neither  small  nor 
large).  The  behavior  theoretically  predicted  for  the 
locus  at  large  p (in  particular  the  rapid  movements 
between  ImX  = n/rev  and  n + l/2/rev,  and  perhaps - 
for  y not  too  large  — between  RHP  and  the  LHP) 
actually  does  show  up  in  the  numerical  calculations  of 
the  stability;  such  behavior  of  a numeric?)  solution 
might  be  questioned  without  the  perturbation  solution  to 
provide  a guide  to  what  to  expect.  It  is  unfortunate  that 
the  boundary  of  the  inb*nbllity  region  for  >/16  of  order 
1 is  first  encountered  at  moderate  (around  p 2.25 
for  small  Kp  ; sou  Fig.  12)  and  so  cannot  be  obtained 
by  perturbation  techniques  (to  the  order  explored  anyway). 
Because  of  the  small  time  constant  in  the  main  solution 
(of  order  t:~ * ) and  the  four  transition  regions  (of  width 
0(p~?  3)  , a numerical  calculation  of  the  roots  for  truly 
large  p would  be  difficult;  the  perturbation  theory 
handles  these  singular  problems  analytically,  and  the 


11 


calculations  that  remain  are  noosin', ular , short,  and 
simple. 

A Comparison  With  Numerical  Results 

Shown  in  Fig.  13  arc  lines  of  constant  ImX  and 
ReX  for  v = 1 and  Kp  = 0 , based  on  numerical 
calculation  of  the  eigenvalues  of  equation  (1)  using 
Floquet  theory  (this  figure  is  taken  from  Fig.  3 of  Ref. 
2).  i gure  13  should  be  compared  with  the  analytic 
results  of  Fig.  3 (small  p)  and  Fig.  12  (large  p)  . 

The  small  p results  are  good  out  to  about  p = 0.5, 
which  is  an  excellent  range,  including  ail  the  operating 
conditions  of  most  helicopters.  The  large  p results 
are  valid  above  p = 2.5  or  so,  which  is  also  a good 
range  for  the  perturbation  theory.  The  first  high  p 
instability  lobf  f Fig.  12  can  be  seen  in  Fig.  13, 
appearing  between  p =2.0  and  2.5 

Application  of  Perturbation  Techniques  to  Helicopter 
Dynamics 

This  section  returns  to  the  question  of  whether 
perturbation  techniques  might  be  profitably  applied  to 
realistic  dynamic  systems  than  the  one  considered  here. 
As  part  of  the  answer,  consider  what  these  techniques 
will  not  do:  obviously  they  cannot  give  results  for  cases 
where  there  is  no  parameter  that  is  either  small  or 
large,  for  example  when  > = 16  and  p = 1 . How- 
ever, the  four  cases  considered  together  cover  a good 
deal  of  the  ranges  of  p and  V , and  with  primarily 
analytical  results.  For  many  helicopters  the  small  p 
case  will  be  quite  satisfactory  alone.  What  the  tech- 
niques can  do  also  includes: 

(a)  Since  they  give  analytic  solutions  they  provide 
more  insight  into  the  problem,  as  well  as  specific  de- 
sign criteria  for  the  system;  this  feature  is  particularly 
important  for  nonlinear  or  time- varying  systems,  which 
have  properties  much  different  from  those  of  constant 
coefficient,  linear  systems. 

(b)  Perturbation  methods  can  find,  and  handle, 
cases  that  are  very  sensitive  to  the  parameters,  or  that 
are  difficult  to  solve  accurately  by  numerical  methods. 

,c)  The  methods  provide  more  insight  into  the 
miner  i nusual  behavior  of  the  solution  of  periodic 
systems,  by  showing  explicitly  how  the  periodic  co- 
efficients modify  the  transient  solutions  and  why  they 
give  the  root  loci  their  characteristic  behavior  in  the 
critical  regions. 

(d)  Finally,  even  if  the  techniques  are  not  used  to 
find  the  complete  solution,  it  only  takes  a little  work  to 
find  out  where  the  problems  are  (e.g. , critical  regions 
and  transition  regions)  and  what  the  order  of  things  is, 
which  information  would  be  of  invaluable  help  in  the 
numerical  analysis  of  a system. 

The  extension  to  more  degrees  of  freedom  or  more 
realistic  aerodynamic  coefficients  would  certainly  make 
the  analysis  more  complicated.  In  general  however  any 
study  — analytic,  computational,  or  experimental  — of  a 


system  becomes  more  complicated  as  the  accuracy  of 
the  modeling  of  the  true  system  increases,  and  perturba- 
tion techniques  are  not  expected  to  be  an  exception  to 
this  rule.  Regardless  of  the  system  being  studied,  the 
position  perturbation  techniques  occupy  between  simple 
linear  analyses  and  complex  nonlinear  numerical  calcula- 
tions makes  them  a very  powerful  tool  for  providing  both 
exact  solutions  and  increasing  understanding  of  problems 
in  rotor  dynamics . 

There  are  many  problems  in  rotor  dynamics  in- 
volving nonlinear  or  periodic  coefficients  to  which  per- 
turbation techniques  might  profitably  be  applied.  Some 
additional  work  might  be  done  with  the  flapping  dynamics, 
considering  for  example  several  blades  responding  to  a 
gust  or  to  shaft  motion.  The  coupled  pitch  /flap  or  flap/ 
lag  dynamics  of  a single  blade  are  other  important 
problems  of  rotor  dynamics  to  which  perturbation  theory 
might  usetully  be  applied.  While  the  perturbation  solu- 
tion for  small  p will  probably  be  of  most  interest, 
there  will  likely  arise  problems  where  other  parameters 
are  also  useful.  As  long  as  a reasonable  model  is 
chosen  for  the  system,  and  as  much  effort  is  given  to  the 
interpretation  of  the  solution  as  to  its  derivation,  per- 
turbation techniques  should  prove  quite  useful  in  provid- 
ing information  about  these  problems,  and  many  others 
in  rotor  dynamics  ard  ac.odynamics. 

T.ds  paper  has  provided  examples  of  the  informa- 
tion about  dynamic  systems  which  may  be  obtained  using 
the  methods  of  perturbation  theory.  The  techniques  have 
proved  very  useful  for  the  problem  studied.  It  should 
not  be  concluded  however  that  the  techniques  used  for 
this  problem  are  all  the’-e  to  perturbation  theory; 
there  are  many  more  metuods  that  have  not  been  touched 
on  here.  Perturbation  theory  is  a powerful,  and  yet  not 
very  sophisticated,  mathematical  technique  which  should 
prove  very  useful  in  analyzing  some  of  the  problems  of 
helicopter  dynamics. 
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Fig.  1.  Foot  loci  for  varying  Y,  baaed  on  the 
small  n results  (to  order  >i2);  v = 1 
and  n - 0 and  0.1. 


Imx 


Fig.  2.  Root  loci  for  varying  n , based  on  the 
small  n results  (to  order  n2);  Y fixed 
for  each  locus  (ReX  = - Y/l6  for 

>i=0) 


Fig.  3.  Lines  of  constant  ImA  and  ReX,  based  on 
the  small  >i  results  (to  order  >i2);  v = 1. 

Kp  = 0;  ImX, ReX, 

circled  values  of  ImX  indicate  areas  in 
which  ImX  is  constant 


Fig.  4.  Lines  of  constant  ImX  and  ReX.  based  on 
the  small  >i  results  (to  order  >j2)  ; » 1 . 

Kp  = 0.1; ImX. ReX. 

circled  values  of  ImX  indicate  areas  in 
which  ImX  is  constant. 


Fig.  5.  Lines  of  constant  ImX  and  ReX,  based  on 
the  small  n results  (to  order  >i2) ; v = I , 

Kp  = -0. 1 ; ImX,—  — ReX, 

circled  values  of  ImX  Indicate  areas  in 
which  ImX  is  constant 


0 , 0 

Fig.  6.  The  ratio  (-Mp/M^),  which  governs  the 

effect  of  K and  K for  small  v 
R P 
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Fig.  7.  The  averages  of  the  aerodynamic  Fig.  8.  Root  loci  for  varying  p and  V,  based  on  j 

coefficients,  which  give  the  roots  the  small  V results  (to  order  V);  v = 1 

for  small  V.  and  Kp  = 1,  the  effect  of  the  critical 

regions  is  not  shown. 


P 

Fig.  9.  The  Kp  vs.  p boundaries  for  stability 
in  the  center  of  the  critical  region,  based 
on  the  small  V results  (to  order  V);  for 
roots  near  ImX  = 1/2,  1,  and  3/2. 


.675 


Fig.  10.  Motion  of  the  root  on  the  real  axis  for  varying 
p , V = «>.  For  large  V there  are  two  real 
roots;  X j is  the  finite  real  root;  the  other 
root  is  at  X = - «o  for  V = *>  . 


NEGATIVE  DAMPING  ON 
RETREATING  SIDE 


Fig  11.  Flap  '•ate  feedback  required  for  negative 
damping  over  part  of  the  rotor  disk; 

- (Mi  - K^M^)  < 0 (enters  into  the 
large  V case). 


H- 

Fig.  12.  Sketch  of  the  characteristic  behavior  of  the 
critical  region  boundaries  and  stability 

boundaries  for  large  p ; 

boundary  of  region  in  which  ImX  is  fixed 

at  n/rev  or  n+l/2/rev; 

boundary  of  region  in  which  the  real  part 
of  one  root  is  positive,  i.e.  unstable 
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Fig.  13.  Lines  of  constant  ImX  and  ReX,  based  on 
numerical  calculation  of  the  eigenvalues 

(from  Ref.  2);  v = 1,  Kp  = 0; 

ImX, ReX;  circled  values  of 

ImX  indicate  areas  in  which  ImX  is 
constant. 
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